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ABSTRACT 

The stellar initial mass function (IMF) describes the distribution in stellar masses 
produced from a burst of star formation. For more than fifty years, the implicit as- 
sumption underpinning most areas of research involving the IMF has been that it is 
universal, regardless of time and environment. We measure the high-mass IMF slope 
for a sample of low-to-moderate redshift galaxies from the Galaxy And Mass Assem- 
bly survey. The large range in luminosities and galaxy masses of the sample permits 
the exploration of underlying IMF dependencies. A strong IMF-star formation rate 
dependency is discovered, which shows that highly star forming galaxies form pro- 
portionally more massive stars (they have IMFs with flatter power-law slopes) than 
galaxies with low star formation rates. This has a significant impact on a wide vari- 
ety of galaxy evolution studies, all of which rely on assumptions about the slope of 
the IMF. Our result is supported by, and provides an explanation for, the results of 
numerous recent explorations suggesting a variation of or evolution in the IMF. 

Key words: galaxies-stellar initial mass function: galaxies-formation and evolution 



1 INTRODUCTION 

The stellar initial mass function (IMF) is an empirical 
power-law relation describing the distribution of stellar 
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masses formed in a single episode of star formation. This 
initial stellar mass distribution has often been assumed to 
be universal. This is perhaps the most fundamental assump- 
tion used in all galaxy formation and evolution studies. The 
IMF is the bridge between the massive stars, measurable 
through tracers such as Ha, ultraviolet, far-infrared and 
radio luminosity, and the low mass stars, which form the 
bulk of the stellar mass in galaxies ( Kennicutt||1998 1. The 
IMF is intimately involved in many aspects of the modeling 
of galaxy evolution. Some examples include models of tur- 
bulent fragmentation and collapse of gas clouds that form 
sub-stellar to super-stellar objects ( [Nakamura fc Umemura| 
2001), the numerical study of supersonic hydrodynamics 



of galaxies and non-traditional approaches, all reporting 
discrepancies between the observations and model predic- 
tions. A number of recent studies now suggest an evolv- 
ing or spatially varying IMF as a "last resort" explana- 



tion to reconcile the observed differences (Hoversten & 



Glazebrookl [2008) [Wilkins, Trentham fc Hopkins) pOOSb 



Wilkins et al.||2008a| [van Dokkum| |2008| [Meurer et al.| 



20091. The assumed cosmological parameters are: Hq 
kms"^ Mpc-\ Q.M = 0.3 and Q,^ = 0.7. 



70 



and magnetohydrodynamics of turbulence (Padoan et al. 
2007|, gradual processes behind building of a galaxy ( |Gib-| 
son fc Matteucci|1997 1 , the reionisation of the intergalactic 
medium at high redshift (2 > 6) ( Chary|2008 1, the relation- 
ship between stellar mass and star formation rate (SFR) 
(Dave 2008 1, evolution in colour and mass-to-light ratio 



2 OBSERVATIONS 

2.1 Galaxy And Mass Assembly survey 

Motivated by these recent failures of the universal IMF as- 
sumption, we have conducted an analysis exploring such 
variations using a sample of galaxies from the Galaxy 



of galaxies Ivan Dokkum 20081, models of heavy element And Mass Assembly (GAM A) survey (Driver et al. 2009 



production, chemical enrichment and evolution in galaxies 
( Calura fc Menci|2009 \ from the death of stars, the fraction 
of stars that form black holes ( Fryer|2003 1 and many other 
evolutionary processes. 

The IMF is often parameterised as one or more power 
laws, describing the number of stars within a given mass 
interval, ^ oc m" with a defining the slope for the 



mass range of interest ( Salpeter|1955 Baldry & Glazebrook 
2003). Other widely adopted functional forms for the IMF 



include a lognormal form for the low mass regime with a 
power-law tail for high masses ( |Chabrier|200l| ). From star 
counts of local resolved stellar populations, the IMF is mea- 
sured to have a high-mass (m > 0.5 Mq) slope of a ~ —2.35 
(Salpeter||1955 Scalo||1986 ), also called the Salpeter slope. 



although variations of this slope are also reported ( Kroupa, 



Tout fc Gilmorel[T9 93': 'Kroupa '2001"Miller fc Scai o|1979 i' 
An "integrated galaxy" IMF (IGIMF) has recently been 
proposed for the interpretation of the galaxy-wide IMF 



properties. Whether the IGIMF is universal (Elmegreen 



2006) or differs from the star cluster IMF (Weidner & 



Kroupa 2005) is again a much debated subject. The de- 



pendence of cluster formation on the SFR of a galaxy is 



argued to give rise to a varying IGIMF ( Weidner fc Kroupa 



2005), although the underlying IMF may still be univer- 



sal. In almost all cases, the high-mass IMF slope is mod- 
elled with a power-law type behaviour and a multi-part 
IMF expression best describes the stellar luminosity func- 
tion in the solar neighbourhood ( Scalo|1986 Kroupa, Tout 
fc Gilmore|1993 ). Therefore, a multi-part power law is used 
to obtain the results presented here. For the purpose of this 
investigation, we primarily use a 2-part power law with a 
Salpeter high-mass (> 0.5 M©) slope and a = —1.3 low- 
mass (0.1 ^ M/Mq ^ 0.5) slope. The use of other popular 



funct i onal forms of the IMF ( [Miller fc Scalo||1979| [Scalo 
1986 Kroupa|[2001 1 does not influence the overall conclu- 
sions of this study. 

Turbulent star-forming gas and clump mass functions 



(Reid fc Wilson 20051 indicate that a Salpeter-like IMF 



slope may be imprinted on the mass distribution of turbu- 
lent structures, and that all stellar clusters formed out of 
these structures therefore inherit a Salpeter-like IMF. The 
concept of a "universal IMF" is, however, being increas- 
ingly scrutinised by recent studies based on large samples 



[Robotham et al.|2009| [Baldry et al.|2009| ). GAMA is a spec 
troscopic survey with multi-wavelength photometric data 
undertaken at the Anglo-Australian Telescope using the 
2dF fibre feed and AAOmega multi-object spectrograph. 
AAOmega provides 5 A resolution sp ectra with complet e 
spectral coverage from 3700-8800 A ( [Sharp et a"L][2006 l. 
GAMA covers three equatorial fields of 48 deg'^ each, with 
two fields reaching a depth of tab < 19.4 magnitude and 
the third extending to r^s < 19.8 magnitude, together 
with Kab < 17.6 magnitude over all three fields. There are 
~ 120 000 galaxies with measured spectra available from 



GAMA observations to date (Driver et al. 2010, submit 



[ted[) . The redshift of each galaxy is determined using RUNZ 



(Saunders, Cannon & Sutherland 2004), a FORTRAN pro 



gram for measuring redshifts from reduced spectra. 



2.2 Data 

The standard strong optical emission lines are measured 
from each curvature corrected and flux calibrated spec- 
trum assuming a single Gaussian approximation and com- 
mon values for redshift and line width (Bauer et al. 2011, 
in prep). Corrections for the underlying stellar absorption, 
dust obscuration and fibre aperture effects, detailed below, 
are applied to these measurements. A full composite line 
and continuum extraction process is ultimately intended 
for the full data set. 

The strength of Ha emission in galaxy spectra is used 
in this investigation to probe the extent of the star for- 
mation in galaxies. The Ha luminosity is used to measure 
the current SFR, as the ionising photons mainly come from 
short-lived massive stars. Our sample is drawn from the 
~ 120 000 spectra available at June 2010, and is comprised 
of 43 668 galaxies with measured emission lines, about 40% 
of all galaxies in the GAMA sample at that time. This sam- 
ple only includes objects with redshift quality flags ^ 3 



(i.e. regarded as a secure redshift, see Driver et al. 2010 



submitted). Furthermore, we exclude all galaxies with Ha 



emission measurements affected by the presence of strong 
sky lines (see Figure[2|, and all galaxies with Ha emission 
below a minimum flux limit of 25 x W^"^ erg/ s/cm^ . This 
flux limit is obtained from examining the spectra of a sam- 
ple of low Ha luminosity galaxies. Increasing this flux limit 
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Figure 1. (a) The distribution of SFR and galaxy stellar masses 
with redshift. The visible gap in the distribution centred on z = 
0.16 shows where the wavelength of the atmospheric O2 band 
(Praunhofer A— line) overlaps with the redshifted wavelength of 
the Ha emission line, leading us to omit these data from our 
analysis. The masses are in units of the solar mass, 1 Mq = 
1.99 X 10^" kg. (b) The definition of the volume limited samples 
used in this study. The galaxies within each magnitude range are 
not affected by the flux limits of the survey, so that each of the 
three samples is complete to a given luminosity. 



to that used by Brough et al. (20101 for a sample of low 
Ha luminosity GAMA galaxies, for example, does not al- 
ter our conclusions. Sloan Digital Sky Survey (SDSS) pho- 



tometry in u,g,r,i,z filters is available for each galaxy (Hill 
et al. 2010 1, fc-corrections to 2 = 0.1 are applied and all 



photometry is corrected for foreground (Milky Way) dust- 
extinction ( Schlegel, Finkbeiner fc Davis|1998 1 . The galaxy 
sample covers a moderate range in redshift (0 < 2 ^ 0.35). 
Galaxies dominated by emission from active galactic nuclei 
(AGN) are excluded from the sample (5334 galaxies) based 
on standard optical emission-line ([NH]/Hq and [OHIJ/H/?) 
diagnostics using the discrimination line of |Kewley et al.] 
(2001 1. In the case of galaxies for which only some of these 
four emission lines are measurable, AGNs can still be ex- 
cluded using the diagnostics log ([NII]/HQf)^0.2 and log 
([OIII]/H/?)^l. This excludes a further 173 galaxies. The 
size of the final sample is 33 657. 

This sample of galaxies spans a large range in stel- 



lar mass (7 sC log (M/M©) ^ 12) and SFR (10" 
Meyr- 



-100 



redshift that permits us to explore the potential IMF de- 
pendencies with respect to different physical properties of 
galaxies. Figure [T]^ a) shows the wide range in SFR sampled 
as a function of redshift, and colour coded to illustrate the 
range in stellar masses. Figurejlj^b) shows the envelope of 
the distribution of absolute r-band magnitude, Mr, with 
redshift for this sample. Outlined within this envelope are 
three independent volume-limited samples which form the 
basis of our subsequent analysis. These are selected to span 
~ 1 magnitude in Mr, centred on the values shown. 



3 DERIVING PHYSICAL QUANTITIES 
3.1 Ha luminosities 

The observed Ha emission must be corrected for the effects 
of stellar absorption within the host galaxy, obscuration by 
dust, and the Umited sampling of each galaxy by the optical 
fibres (aperture effects) used in multi-object spectroscopy. 



3.1.1 Stellar absorption correction and Balmer 
decrements 

A simple constant correction for stellar absorption in 
Balmer emission line equivalent widths (EWs) (i.e. Ha and 
H/3 EWs) is used for this investigation. The assumed com- 
mon EW correction (EWc) for the stellar absor ption in the 
GAMA d ata is 1.3 A. Based on previous work dHopkins et 
al. 20031, a correction of at most 1.3 A is sufficient, pro- 



vided the assumption is restricted to studies examining the 
gross characteristics of a large sample of sources, which is 
the case in this investigation. We tested a range of EWc 
values between 0.7 and 1.3 A, and the results did not vary 
measurably. Only Ha EWs smaller than log (Ha EW)< 0.9 
show a difference in EW of more than 5%. 

As the stellar absorption may in general be a lumi- 
nosity dependent quantity, the required correction could be 
higher for high star formation rate sources. In order to test 
this aspect, a unique EWc for each galaxy is assigned based 
on an assumed linear relationship between luminosity and 
stellar absorption correction. The results indicate that a lu- 
minosity dependent stellar absorption correction does not 
significantly affect the calculated intrinsic luminosity of the 
source. Less than 10% of the sample showed any noticeable 
effect. For systems with log (Ha EW)> 2, the assumption 
of a luminosity dependent absorption correction increases 
the inferred EWs for these extreme systems, thereby fur- 
ther enhancing the trend with SFR presented in this paper. 
Therefore, the assumption of a fixed EWc = 1.3 A should 
not significantly affect the trends evident in the results, or 
our conclusions. 



The Balmer Decrement (BD) is defined as the ratio 
of stellar absorption corrected Ha to H/3 fiuxes {BD = 
Shci/ Shp), where Sna for example is ( Hopkins et al.|2003 l. 



{HaEW + EWc) 
HaEW ' 



(1) 



^). It is this large range in SFR, stellar mass and 



where Fhu is the measured emission line flux. 

A small fraction of galaxies have Balmer Decrements 
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Figure 2. (a) The relationship between Balmer Decrements and 
aperture corrected luminosities for the galaxies with accurate Ho 
and H/3 EW measurements. A linear relation is fitted to the data 
to determine the Balmer Decrements for the galaxies that have 
measured H« EWs but no H/3 EWs or have H/3 EWs that are 
affected by sky absorption bands, (b) The dust obscuration as 
indicated by the Balmer Decrement is largest at high z. This is 
because high SFR objects also tend to have higher obscurations 
( [Hopkins et al.||2001| and the objects at higher z tend to have 
higher SFRs, a result of both galaxy evolution and a predomi- 
nance of massive, high— SFR galaxies at high-z, due to the flux 
limit of the survey. 



(BD) less than the Case B recombination value BD= 2.86 
(Figure |2|. The Balmer Decrement is an obscuration sensi- 
tive parameter, and its departure from the Case B recombi- 
nation value of 2.86 is an indication of the dust attenuation 
along the line of sight. BD< 2.86 can result from an in- 
trinsically low reddening combined with uncertainty in the 
stellar absorption, but also from errors in the line flux cali- 
bration and measurements ( Kewley et al.|2006 1 . Out of all 
the galaxies with measured Balmer Decrements, ~ 6% have 
BD< 2.86. All these Balmer Decrements are set to 2.86 for 
the purpose of this investigation. 

Not all galaxies have both Ha and H/3 measurements. 
For the galaxies with only Ho? measurements, the relation 
between aperture corrected luminosity and Balmer Decre- 
ment is used to determine the Balmer Decrements. The em- 
pirical form of the relationship between aperture corrected 
luminosity and Balmer Decrement is shown in Figure [2] 




0.35 



Figure 3. Aperture correction factor as a function of redshift. 
The required aperture correction is largest at low z because low 
2 objects are more likely to be larger in angular size than the 
aperture of the spectroscopic fibre used for the observations. 



The form is 

BD = 0.811og(L/f<,,ACor 



2.86, 



22.041, log{LHc.,ACor) > 30.74 

log{LHc,ACor) ^ 30.74 (2) 



where \og{LHa,ACor) denotes the aperture corrected Ha lu- 
minosity (see Eqn.|3|. 



3.1.2 Aperture correction 

Aperture effects arise from the physical limitation imposed 
by the diameter of the spectroscopic fibre used in the ob- 
servations. For nearby sources, this means that the fibre 
only captures part of the light from the object, which is 
naturally a problem for sources larger in size than the fi- 
bre diameter projected on the sky. An aperture correction 
is required to account for the missing flux, in order to get 
an estimate of the true star formation rate. Following the 



approach used for SDSS spectra by Hopkins et al. (20031 



we implemented an aperture correction for the Ha lumi- 
nosities of the GAMA galaxies. This approach uses the ab- 
solute 1 — band magnitude to approximate the continuum at 
the wavelength of Ha, thereby accounting for the Ha lu- 
minosity of the whole galaxy (Eq.[3|. Figure[3] shows the 
relation between the applied aperture correction and z and 
the required correction is typically a factor of 2 — 4. 

The aperture corrected Ha luminosity {LHa,ACor) for 
the whole galaxy is a function of three parameters. Ha EW 
(EWHa.obs), absolute r-band luminosity [Mr) and redshift 
(z). The form of the aperture corrected luminosity before 
applying any obscuration correction is I Hopkins et al.|20^ I 



-^HctjApCor 



(Wh.,...+Wc)x10-°-^(*^'^-3^-^°) 

3 X 10^** 
^ (6564.61(1 + 2))2' 



(3) 



The aperture corrections are based on the absolute 
magnitude, M,-, of each galaxy as an estimate of contin- 
uum luminosity, thereby recovering a Ha luminosity for 
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Table 1. The SFR calibration factors for different IMFs 



high- 


-mass slope 


calibration factor (W) 


a 


= -2.00 


4.08 X 1034 


a 


= -2.35 


1.27 X 10^4 


a 


= -2.50 


0.65 X 10^'' 



the whole galaxy. This method of applying aperture correc- 
tions to the luminosities, described in Hopkins ct al. ( 20d3| , 
yields similar results to the more complex colour gradient 
based aperture corrections described in |Brinchmann et al.| 
(20041. We use these aperture corrected values throughout 



this analysis. We have, in addition, tested the efTect on our 
results in the case of using SFRs estimated only from the 
detected Ha emission through the fibre (no aperture correc- 
tions). Even in this case, we find the same qualitative con- 
clusions regarding the SFR-dependence of IMF slope. We 
conclude that the aperture corrections applied here are not 
introducing any significant bias, nor are they erroneously 
giving rise to our results. 

The aperture, obscuration and stellar absorption cor- 
rected luminosity for the whole galaxy is given as 



3 X 10^** 



2.86 



(6564.61(1 + 2))2 V 2.86 7 ' ^ ^ 

where Fhc/ Fnfi denotes the Balmer decrement. Figure|2|b) 
explores the increase in Balmer decrement with respect to 
increasing z. The exponent of the Balmer Decrement in 
Eq.|4]is equal to k{\Ha)/[k{\Hp) — ^(Ai/a)], where extinc- 
tion at a given A, kiX), is determined from the |Cardelli,| 
Clayton & Mathis ( 1989 1 galactic dust obscuration curve. 



The Hq SFR can be determined from Kennicutt ( 1998 1 



SFRhc 



1.27 X 1034' 



(5) 



The derivation of the Ha SFR calibration requires the 
assumption of an IMF. The above calibration factor has 
been derived assuming a Salpeter IMF. Table[l] shows the 
effect on the SFR calibrator if a different IMF is assumed. 
For a given luminosity, a calibration based on a fiatter IMF 
would indicate a lower SFR than the Salpeter IMF based 
calibration. 

If an IMF dependent SFR calibration is used in the 
derivation of SFRs for GAMA galaxies, that would reduce 
the range in SFR shown in Figure 1(a). This reduction in 
range would not affect our main conclusion of an IMF- 
SFR relationship, because the SFR calculated for the sam- 
ple would still vary monotonically (as the scaling is linear, 
the ordering of the SFRs is not affected). 



3.2 The determination of stellar masses 

The stellar masses used in this investigation are derived 
based on the observed tight relation between [g — i) colour 



and the mass-luminosity (M/L) relation (Taylor et al. 
2011; in prep) using the Bruzual & Chariot (20031 



2010 



models and assuming a Chabrier (20031 IMF. This method 



of calculating the masses yields results consistent with other 
established techniques (e.g., Baldry et al.||2006 l. 



log(M,) = -0.68 + 0.73(g - i) - QA(Mi - Af^,©), 



(6) 



where M* is the mass of the galaxy, g and i band colours 
are fc-corrected to z = 0, Mi is the absolute magnitude of 
the galaxy in the i-band and Mi^Q — 4.58, the absolute 
magnitude of the Sun in the i-band. 



4 OBSCURATION CORRECTIONS 

Understanding and interpreting the physical and chemical 
properties of galaxies depends in part on how accurately the 
data are corrected for stellar absorption and dust obscura- 
tion to recover the intrinsic fiuxes. Hq can be heavily atten- 
uated by dust. High SFR galaxies are subjected to greater 



dust obscuration than lower luminosity objects (Hopkins 
et al.|2001 , 2003] [Afonso et aT][2003l [Perez-Gonzal ez et al 



2003 1. Obscuration corrections are especially critical in this 
analysis as our primary aim is to compare the observed Ha 
EW and g — r [or g — i) colours with PEGASE generated 
synthetic spectra for different input IMFs, assuming no ex- 
tinction. 

The reliability of the applied dust correction depends 
on the adopted dust obscuration models. We explore several 



popular empirical dust models ( Cardelli, Clayton & Mathis 



1989 Calzetti 2001 Fischera & Dopita 20051 along with 



radiative transfer model predictions of the effects of dust 
extinction (Popcscu ct al.|,2000 Tuffs et al.||2004 Popescu 



eraL][20Tl| ). 

The differential reddening between the stellar contin- 



uum and gas (Calzetti 20011 must be addressed in deriv- 



ing the intrinsic fiuxes. The difference in attenuation be- 
tween gas and continuum is generally assumed to be ~ 2 
(Wijesinghe et al.|[2010| |Meurer et al 



2009 Hoversten & 



Glazebrook 2008 C alzetti||2001[ ). [Hoversten fc Glazebrook 



( 2008 1 describe the effect, on Ha EWs and colours, of vary- 
ing the differential reddening factor between the continuum 
and gas. We tested the impact of this assumption on our 
results. For all the subsequent analysis in this paper we use 
obscuration corrected colours derived through the applica- 
tion of an obscuration curve together with this factor of ~ 2 
(as detailed below). These measurements were compared 
against colours derived by Taylor et al. (in prep), from full 
SED modelling of the GAMA photometry, with indepen- 
dent dust corrections. The results are consistent, with an 
RMS scatter of ~ 0.17 mag, with no systematic deviation, 
as might be expected if the factor of ~ 2 between gas and 
continuum obscuration were significantly in error. There is, 
moreover, no systematic offset in these two approaches as 
a function of SFR, specific SFR, mass or redshift. We con- 
clude that our assumption of this commonly used factor is 
justified, and unlikely to introduce any systematic error in 
our result. 

The obscuration corrected Ha EW is given as 



EWh 



OL.int — 



{EWHc.obs + 1.3) 
(1 + 2) 

y^lQOA[k(\Hc.)E(B-V)gas-k(\Ha.)E(B~V),]j^^ 
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where k(A) gives the extinction at wavelength A, and the 
colour excess of gas (i.e., emission lines) is 



E{B - V)gas 



log (if /2.86) 



(8) 



OA[k{XHp)-k{\Hc,)]' 
The colour excess of the continuum ( Calzetti|1997 1 is 

E{B~V),=O.UE{B-V)sas. (9) 

The effect of the spectroscopic fibre sampling only the 
central regions, for those galaxies largest on the sky, may 
be to limit the detection of the lowest EW systems in low 
SFR galaxies. Accounting for this effect would not change 
our results; indeed, if anything, such a correction would act 
to enhance the trend investigated here. 

4.1 Applying obscuration corrections 

This section describes the different obscuration curves used 
in this analysis, and the methods of applying obscuration 
corrections to the data. 

We have used a combination of the |Calzetti| ( |2001[ ) and 



Cardelli, Clayton & Mathis| ( 1989 1 obscuration curves and 



the Fischera & Dopita ( 2005 I curve as described by Wi- 
[jesinghe et al.| ( |2010[ ) to determine the necessary correc- 
tions. These extinction curves and the application of dust 
corrections are described below. 



4-1.1 Calzetti (20011 obscuration law 



This dust extinction curve is appropriate for continuum at- 
tenuation corrections as this curve is derived from spatially 
integrated colours of the entire stellar population in a sam- 
ple of starburst galaxies. Embedded within the analytical 
form of this curve are dust geometry and composition, and 
it mostly describes dust absorption, since the effects due 
to scattering are averaged out. Because the entire stellar 
population within the galaxies was observed, emission lost 
through dust scattering out of the line of sight is averaged 
out by the scattering into the line of sight. The form of the 
curve is 



fc(A) = 2.656 - 1.857 + 



1.040 
A 



+ 4.05, 



for 0.63/im ^ A ^ 2.2/im 
,/ 1.509 



2.656 



2.156 



0.198 0.011 



A 



A2 



A3 



for 0.12^im ^ A < 0.63/im. 



4.05, 
(10) 



The second term of the exponent in Eq.[7| which is 
related to the continuum luminosity, and the corrections to 
the colours are derived using this curve. 



4.1.2 Cardelli, Clayton &: Mathis (1989) 



obscuration law 

This Galactic dust obscuration curve is derived from obser- 
vations of the UV extinction of stars, as well as using various 



other sources for optical and NIR data. The young stellar 
populations in massive star forming regions responsible for 
UV radiation are also responsible for nebular emission lines 
and this curve is applicable to both diffuse and dense stellar 
regions. This curve accurately describes the dust effects on 
emission lines and has the form 



fe(A) = a{x) + 



Rv 



(11) 



where Rv is the ratio of total to selective extinction and is a 
constant for a given extinction curve. The value of i?„ — 3.1 
( |Calzetti|2001 1, which is found to well describe the redden- 
ing of the ionised gas in star forming galaxies, is used in 
this analysis. 

The functional forms of a{x) and 6(2;) in Eq. |ll[ with 
X — 1/A, are a power law in the infrared regime and a 
polynomial in the optical/NIR regime, in units of fim~^: 

Infrared; 0.3/im~^ ^ 2; ^ l.l/xm~^; 

a{x) = 0.574A^-''^ 
b{x) = -0.527A^''\ 

Optical/NIR: 1.1/tm"^ a; s: 3.3/im"^ and y = (a; - 1.82); 

a{x) = l + 0.17699y-0.50447y^- 0.02427/ + 0.72085/ 

+0.01979/ - 0.7753/ + 0.32999/; 
b{x) = 1.41338y + 2.28305y^ + 1.07233/ - 5.38434/ 



-0.62251/ + 5.30260/ - 2.09002/. 



(12) 



The first term of the exponent in Eq.[7] related to the 
line luminosity is based on this curve. The intrinsic Hq EW 
is the ratio of the corrected Hq line to continuum luminosi- 
ties. 



4.1.3 [Fischera &: Dopita] ( [2005] ) obscuration law 



A recent study by Wijesinghe et al. (20101, looking at dust 



obscuration in galaxies using GAMA data, tested a num- 
ber of common obscuration curves, including t he |Calzetti| 



( 2001[|1997| , and ICardelU, Clayton fc M athis (1989) dust 
curves. They found that a Fischera & Dopita (2005) obscu- 



ration curve with R„ = 4.5 and the 2200Abump removed 
gives an excellent agreement between far ultraviolet, near 
ultraviolet. Ha and [OH] derived star formation rate indi- 
cators. 

In this case, both terms of the exponent in Eq.[7] and 
the corrections to the colours are determined using the|Fis-| 



chera & Dopita ( 2005 1 curve. 



4.2 



Popescu et al. (2000 2011 1 and Tuffs et al 



(2004) radiative transfer models 



In addition to the dust corrections based on the above dust 
obscuration curves, the effects of dust attenuation on Ha 
EW and g — r parameters can be determined using radia- 
tive transfer models fPopescu et al.|2000[ [Tuffs et al.|2004 



Popescu et al. 2011 1, where the attenuation of star light 



from disk galaxies with different dust geometries of different 



© 2002 RAS, MNRAS 000,[T|fT8l 



IMF-SFR Relationship 7 



stellar ages constrained by UV/optical to FIR/submm spec- 
tral energy distributions are considered. The model predic- 
tions are based on the opacity of the diffuse dust compo- 
nent, given as face-on B-band optical depth (rb), inclina- 
tion of a galaxy and a dumpiness factor (F) describing the 
local absorption of UV light from massive stars due to the 
presence of massive star forming regions. 

These various approaches to obscuration corrections 
are detailed as vectors in Figure[4] showing the effect of 
different obscuration curves, or models, on the data, for an 
assumed Balmer decrement of BD = 4. As shown in more 
detail below, the different approaches to dust correction do 
not change our qualitative conclusions. 



5 EVIDENCE FOR A NON-UNIVERSAL IMF 

The three model evolutionary tracks shown as black lines 
in Fig ure [5| are reproduced from |Hoversten fc Glaze-| 
brook ( 2008) . These evolutionary tracks ar e generated us- 



5.1 The effect of dust extinction 



ing the population synthesis code PEGASE ( Fioc & Rocca- 



Volmerange 19971, under the assumption of no extinction 



and an exponentially declining star formation history with 
an e-folding time of 1.1 Gyr. The model evolutionary tracks 
denoted by the red lines are generated by combining |Maras-| 



ton 



(20051 and PEGASE models. Based on a similar anal- 



ysis of Ha EW and g — r colour, |Hoversten fc Glazebrook] 
( 2008] ) suggest a possible systematic variation in the IMF 
slope, in which faint galaxies prefer steep IMFs. Our results, 
(Figure [5| shown for three sub-samples based on SFR, are 
consistent with those of Hoversten & Glazebrook (20081, 



and are not sensitive to the choice of population synthesis 
models. The low luminosity systems, which also have low 
SFRs (Figure |5| a)), lie below the central model track rep- 
resentative of a Salpeter IMF (a = —2.35) and towards the 
bottom track with a = —3. In contrast, those with high 
SFRs lie towards the top model track with a = —2. 

A clear variation in IMF slope with the SFR of the host 
galaxy is evident in Figure [5j where the high star formation 
rate systems are characterised by a flatter IMF. In order to 
quantify this effect, the same analysis is performed using 
three independent volume limited samples, with absolute 
r-band magnitude ranges centered on M,. = —21,20.5 and 
— 19.5 shown in Figure [ijb), where each sample is com- 
plete to a given r-band luminosity. This avoids the bias 
against lower SFR systems at higher redshifts, imposed by 
our optical/near-infrared magnitude and Ha flux-limited 
selection. Each sample is further divided into eight sub- 
samples based on SFR (Figure [6|. The clear progression 
towards a top-heavy, or flatter, IMF with increasing SFR 
is evident in all three independent volume limited samples. 

The general trend measured here is that low SFR 
galaxies populate the lower right of the Ha EW and colour 
plane, being characterised by Salpeter, or steeper, IMF 
slopes. With increasing SFR, the galaxy population moves 
upwards to the left from a steep to a flat IMF track, imply- 
ing a SFR dependence of the high-mass IMF slope. Given 
the low redshift range (0 < z ^ 0.35) of the sample, the 
observed IMF-SFR effect cannot be a result of merging 
systems since the merger rate is very low, ~ 2% (De Pro- 



pris et al.|2007 Lotz et al.|2008 ) at these low redshifts, and 
the high-SFR systems are not dominated by mergers. 



FigurejT] compares the uncorrected results with dust ob- 
scuration corrected using the [Cardelli, Clayton fc Mathis] 
( 1989| and Calzetti ( 2001 1 obscuration curve combination 



and the 
given by 



Fischera & Dopita (,2005) obscuration curve as 



Wijesinghe et al.| \2Q1Q\ 



The effect of dust from each of the prescriptions inves- 
tigated is to move data points parallel to the model evolu- 
tionary tracks. The progression of data orthogonal to the 
tracks, as SFR varies, is unlikely to be a consequence of er- 
roneous dust obscuration corrections. Finally, increasing rj 
does not affect Ha EW except at high inclinations, where 
Ha is attenuated more than the r-band due to the low scale 
height. Therefore, if there is an increase in disk opacity (tj,) 
as a function of SFR, this effect would not cause a sys- 
tematic shift of data points orthogonal to the model tracks, 
mimicking the trend with SFR presented in this paper. Any 
such effect due to the presence of dust causes the raw data 
points to move downwards parallel to the model tracks. 



5.2 Addressing the systematics 

Here we explore how the modification of the other free pa- 
rameters of the population synthesis models, as well as the 
introduction of additional models, affects our results. 



5.2.1 Effects of modifying the free parameters 

The three model evolutionary tracks plotted in previous fig- 
ures and shown as dashed lines in Figure[8]are generated as- 
suming an exponentially decreasing star formation history 
with e-folding time of 1.1 Gyr, Z — 0.02 (solar metallicity) 
and an upper stellar mass limit of 120 M© for different in- 
put IMFs. Figure|8] explores the effects of varying these free 
parameters in generating model evolutionary tracks and the 
use of other models. 

The modification of the free parameters of PEGASE 
model does change the position of the evolutionary tracks. 
These changes however, largely move the tracks towards the 
bottom of the plots. Although this can potentially explain 
data in lower-left of the diagram (i.e. low SFR galaxies) 
without resorting to evolving IMFs, in none of the cases 
can it explain the data in the upper right of the diagram, 
as shown in Figure|8|a-c). In other words an evolutionary 
track with an input flat IMF is required to describe the 
highest SFR sources, while the positions of low SFR galax- 
ies can be explained either by using evolutionary tracks 
with input Salpeter-like or steep IMFs and r = 1.1 Gyr or 
with a top-heavy IMF with varying e-folding times. 



5.2.2 



population synthesis models 



Figure[8|d) shows the positions of the Maraston ( 2005 1 
model evolutionary tracks and the highest SFR objects in 
the sample with respect to the PEGASE tracks. The Maras- 
ton tracks are generated based on the same input criteria 
used to generate the three PEGASE tracks. In this part 
of the analysis we investigate the effect of the inclusion of 
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Figure 4. Distribution of all GAMA galaxies up to z=0.355, after dust corrections as given in the text. All data and the model tracks 
are fc-corrected to z=0.1. The colour contours indicate the data density and the three solid lines indicate the three different evolutionary 
paths a galaxy would take if all star clusters within that galaxy have an IMF with a slope of a = —3 (bottom track), a = —2.35 (middle 
track) or a = — 2 (top track). These model tracks are generated using PEGASE. The arrows depict the dust vectors. The red arrows 
represent radiative transfer model predictions calculated using the model of |Popescu et ah] | |2000[ |2011| and [Tuffs et al.] |[2004) and 
from left to right correspond to tj, = 8, 4, 1, all assuming a median galaxy inclination of 60° and F = 0.35. The rest of the vectors show 
the movement of data points for different dust extinction curves and for a Balmer Decrement of 4. Blue: The dust vector calculated 
using the Calzctti (1997) curve for the continuum corrections and CardcUi, Clayton & Mathis (1989) curve for emission line corrections. 
Green: The dust vector corresponding to corrections calculated using |Fischera fc Dopita, (,2005) curve as modified by [Wijesinghe et al.| 
( |2010[ l . Black: The dust vector corresponding to the |Cal2etti| ( |2001[ l and [Cardelli, Glayton &: Mathis] | |1989[ l curves for the continuum 
and emission corrections respectively. 




g-r 

Figure 5. The sample of GAMA galaxies divided into three sub-samples based on SFRs. (a) <SFR (Meyr"^) < 3, (b) 3 sgSFR 
(M0yr~^) < 13 and (c) SFR (MQyr~^)^ 13. The two sets (black and red) of three solid lines indicate the three different evolutionary 
paths a galaxy would take in Ha EW and g — r colour if all star clusters within that galaxy have an IMF with a slope of a = —3 
(bottom track), a = —2.35 (middle track) or q = —2 (top track). The black lines are the evolutionary paths predicted by PEGASE 
Fioc &: Rocca-Volmerange|1997| and red lines are paths predicted by |Maraston| ( |2"005| models. The age increases along the tracks from 
100 Myr (top left) to 13 Gyr (bottom right). Coloured contours are drawn based on data densities of each sub-sample. The ranges in 
data densities are indicated alongside the colour bars of each plot. A representative uncertainty on individual measurements is indicated 
by the error bars in the bottom left of (a). A variation with SFR is apparent across the three panels, with high star forming sources 
evidently preferring a flatter IMF. 
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Figure 6. Three independent volume limited samples {Mr, z) = -21,0.29 (a), {Mr, z) = -20.5,0.19 (b), and {Mr,z) = -19.5,0.13 (c) 
divided into eight SFR sub-samples. A representative uncertainty on individual measurements is indicated by the error bars in each 
top— left panel. The IMF-SFR relationship is evident in each independent volume limited sample. The coloured contours are based on 
the data densities and the ranges of these densities are indicated alongside each colour bar. 
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Figure 7. (Mr,z) = —21,0.29 sample divided into eiglit SFR sub— samples, (a) Data corrected for obscuration using [Cardelli, Clayton] 
fc MathisI l|1989|l fo r emission line correctio ns and|Calzetti| ||2001|l fo r continuum corrections, (b) data corrected for obscuration using 
Fischera &: Dopita| l |2005[ l as prescribed by [Wijesinghe et al.| ( |2010"] l and (c) data without any obscuration correction applied. Model 
tracks as in previous figures. 
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Figure 8. (a) The effect of varying t values on the model evolutionary track a = —2.35 with r = 1.1 Gyr. r = 0.1, 0.5, 0.8, 1, 1.5, 2.5, 13 
Gyr e-folding times are shown as solid lines. The e-folding time increases from left to right and all solid tracks have input IMF slope of 
a = —2.35. (b) The effect of metallicity on the model tracks. The two tracks with solid lines have a = —2.35, the same as the middle 
dashed lined track, which has Z = 0.02, i.e. solar metallicity. The solid line extending above the dashed line at high g — r and low Ho 
EW has Z = 0.05 and the other has Z = 0.005. (c) The variation in the universal IMF track, i.e. a = -2.35 (0.5 < M/Mq < 120), with 
respect to two different assumed stellar upper mass limits. The stellar upper mass limits of the three tracks are 120 (dashed track), 
100, 90 Mq (solid lines), from top to bottom, (d) Solid lines denote Maraston model tracks. The dashed tracks in all panels are the 
three main tracks shown in previous figures. The time spans of all of the tracks are 100 Myr to approximately 13 Gyr, running from 
top left to bottom right. Also shown in each panel is the positions of the high SFR (> 13 Mq yr~^) galaxies in the GAMA sample. 



the Thermally Pulsating-Asymptotic Giant Branch (TP- 
AGB) phase on the model evolutionary tracks. Note that 
the TP-AGE phase is included in the PEGASE models; 
however, it may be the case that PEGASE models may 
not be adequate at describing real stellar populations with 
TP-AGB stars ( |Maraston|2005| . 



The current version of the Maraston ( 2005 1 models 



only provides continuum fluxes and thus colours for a given 
IMF, metallicity and star formation history. The Ha EWs 
are calculated by combining the Ha line flux from PEGASE 
with the continuum flux from the Maraston models at each 
common time step. The presence of TP-AGB stars has a 
signiflcant effect on the colour of a galaxy; however their 
contribution to Ha emission is expected to be negligible. 
The turn-over of the Maraston tracks evident in Figure|8] 
occurs at around the age (0.2 ^ t/Gyr ^ 2), when the 
TP-AGB contributions become important. The difference 
in model tracks at this age arises due to the higher contin- 
uum fluxes given by the Maraston model in comparison to 
the PEGASE model. Despite this difference at early times, 
it is clear that, even using the Maraston models, the higher 
SFR systems favour flatter IMF slopes. 



We further te sted the STARBURST99 models ( |Lei^ 
therer et al.||199"9 l to explore whether there is a variation 
in Ha emission predicted by different models. We found 
that the use of STARBURST99 or PEGASE gave essen- 
tially identical results for Ha emission. We conclude that 
the choice of population synthesis models does not signifi- 
cantly alter our main results. 



distributions. The GAMA galaxy distribution is shown in 
Figure[9|a), and the majority of the sample indicates a pref- 
erence for a flatter IMF. This is a consequence of the dif- 
ferent redshift distributions of the two galaxy surveys, as 
demonstrated in Figure |9];b,c). GAMA has (z) « 0.2 within 
this range, while SDSS has {z) x 0.08. The GAMA sam- 
ple is dominated by relatively high z galaxies, with higher 
SFRs, and which we have demonstrated are those that 
favour a flatter IMF slope. 



6 STAR FORMATION BURSTS 

A sudden burst of star formation on top of an otherwise ex- 
ponentially declining star formation history would give rise 
to a large EW (i.e. increased SFR) and make the galaxy ap- 
pear blue for a short period of time. The effect of a burst is 
therefore to push the EWs of galaxies with a Salpeter IMF 
to high Ha EW and low g — r, potentially leading to the 
erroneous inference of a flatter IMF for such a galaxy if only 
its position in the Ha EW and g — r plane is considered. 
This was explored in some detail by [Hoversten fc Glaze-| 
brook] ( |2008] ), who argue bursts are unlikely to explain the 
variation in IMF, as the necessary bursts would have to be 
unrealistically coordinated in time to produce the observed 
galaxy colours. Here we describe several methods that we 
employed to rule out bursts as the possible source of the 
observed IMF-SFR dependency. 



5.3 SDSS vs GAMA 



A recent study by Hoversten & Glazebrook (20081, using a 
sample of ~ 130 603 SDSS galaxies spanning < z ^ 0.25, 
found that the majority of the galaxies in their sample, 
which is dominated by low-luminosity systems, prefer steep 
IMFs. The same redshift limits are used with the GAMA 
sample in order to compare the GAMA and SDSS galaxy 



6.1 Mass— doubling times 

The mass-doubling times (td) of galaxies provide a method 
of isolating those galaxies undergoing a star-bursting phase. 
The mass-doubling time, used to calculate the id values 
for the highest SFR sub-group of galaxies, is deflned as 



ti = M./(0.5 X SFR) ( [Noeske et aL][2007| ), and provides 
a time scale within which the current SFR would produce 
the observed stellar mass (M«). Folded within the constant 
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Figure 9. (a) The distribution of the GAMA galaxy sample spanning < ^ ^ 0.25. The model tracks from the top to bottom have 
the following IMF slopes, a = —2, —2.35, —3. This figure is similar to Figure 1 of |Hoversten fc Glazebrook| l |2008l l; the scale and model 
IMF tracks of both figures are the same, (b) The distribution of the redshift of the SDSS main sample and the GAMA sample for 
< 2 0.25. (c) The GAMA sample colour coded according to z. This sample covers the range < z ^ 0.25. It is clear that the 
reason the peak in the data density in Top-panel lies higher than in the corresponding diagram from [Hoversten fc Glazebrook| | |2008[ l 
is because of the higher redshift range sampled. 
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Figure 11. (a) The highest volume limited sample Mr = —21 divided into sub— groups of 'specific' SFR (S denotes the specific SFR 
in units of yr^^). Specific SFR increases from left— to-right, top-to-bottom, and (b) summaries this trend clearly. 



term (1 — 7? = 0.5) is the IMF dependent gas recycling fac- 
tor, which determines the fraction of gas recycled into the 
interstellar medium. Galaxies with td significantly shorter 
(td < 0.1 — 1 Gyr) than the adopted ages of galaxies are 
potential bursts. A galaxy experiencing a burst of star for- 
mation has a relatively high specific SFR (i.e. SFR per stel- 
lar mass) and hence corresponds to a low t^. The derived 
td values for the highest SFR sub-group presented in Fig- 
ure 10 a) show that the majority of high-SFR galaxies have 



relatively large td values (td > 1 Gyr), indicating that they 
are not currently in a starburst mode. In addition, both t^ 
and mass vary smoothly along the PEGASE model tracks 
from top-left to bottom-right, without the vertical gradi- 
ents in colour that would be expected from starbursts. The 
IMF dependence in the estimation of td will be small. Since 
td oc M/ SFR, the dependencies of mass and SFR on the 
IMF largely cancel (although not entirely, as these depen- 
dencies are not identical). The dependency of R on the IMF 
is likely to have a small effect on the td values, as R is in- 
versely related to the slope of the IMF (e.g. R « 0.28 for 
a « —2.35 and R ~ 0.56 for a ~ —2. 15; [Hopkins fc Beacoin] 
2006 1, such that td increases if a fiatter IMF is assumed and 



vice versa. 

According to the staged galaxy formation scenario of 
[Noeske et "aT] ( |2007[ ), the high specific SFRs (low td) of 
most low mass galaxies (MsJ Mq) are not indicative of 
evolved galaxies experiencing a starburst. In fact an initial 
burst followed by gradual decline seems to be the favoured 
mode of star formation. However, this initial dominant 
burst of star formation is pushed towards later redshifts for 
less massive galaxies. This further supports our argument 
that the highest star formation rate objects of our sample 
are not starbursts but quiescently evolving galaxies inher- 
ently preferring a flatter IMF in comparison to low SFR 
systems. We reiterate that the smooth IMF-SFR variation 
is evident in all three independent volume limited samples. 
The observed trend is thus strong evidence for an IMF-SFR 
relation. 

We next measure the dependencies on specific star for- 



mation rate and star formation rate surface density, to allow 
an exploration of the most likely underlying dependency of 
the IMF. 



6.2 Specific star formation rates 

Specific SFR is calculated as the SFR per unit stellar mass. 
Figure fTT] shows the variation with specific SFR for the 
highest redshift volume-limited sample. This demonstrates 
a very similar result to that found above for absolute SFR, 
in the sense that systems with higher specific SFR also pre- 
fer flatter IMF slopes. Note that even the highest specific 
SFR systems here still have mass doubling times of the or- 
der 1 Gyr, a consequence of their high average masses (see 
also Figure 10 1. The smooth decline along the model evo- 



lutionary tracks with decreasing specific SFR, evident in 
Figure fTT] is consistent with smooth star formation histo- 
ries. Similar results are seen in the other two volume limited 
samples. 



6.3 Star formation rate surface density 

Figure [T2] explores the trends with respect to SFR per unit 
surface area in M0yr~^kpc~^. The Petrosian radius, which 
is derived from the surface brightness profile of the galaxy, 
provides a measure of the angular size of the galaxy. As- 
suming the Petrosian radius to be the radius of a circle, 
the surface area projected on the sky can be calculated, 
from which the SFR per surface area can be determined. 
However, galaxies are not typically circular, and there is a 
range of morphological types, with ellipticals, spirals and 
irregulars all present in the local galaxy population. The 
surface areas projected on the sky by these types are best 
described using ellipses. The effect of assuming a circular 
form is to reduce the SFR per area. This results in mov- 
ing objects between the SFR surface density bins shown 
in Figure |12| Given the large ranges of SFR per surface 
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Figure 12. The highest volume limited sample Mr = —21 di- 
vided into sub-groups of star formation rate surface density 
(i.e. SFR per unit area, in units of M0yr~^kpc~^). The SFR 
surface density increases from left to right, top to bottom. 



area of the eight bins, the interchange of objects between 
bins is unhkely to be significant, and is dominated by ob- 
jects within the scatter outside the lowest density contour. 
The central contours showing the highest density regions 
are not affected by a low interchange of objects. Hence, we 
can categorically say that all objects in the highest SFR 
per area bin have high SFR surface density and certainly 
prefer a shallow IMF slope. Furthermore, the SFR-IMF de- 
pendency shown is not primarily due to galaxy sizes, e.g., 
massive galaxies having higher SFR due to the SF processes 
being distributed over a larger area than the less massive 
systems. The other volume limited samples follow similar 
trends as shown in Figure[T2]for the highest redshift volume 
limited sample. 

These results demonstrate a qualitatively similar de- 
pendence of the inferred IMF slope on SFR, SSFR and SFR 
surface density. The following section quantifies these de- 
pendencies, allowing us to identify which is likely to be the 
more fundamental. 



7 THE FUNDAMENTAL IMF DEPENDENCY 

Here we attempt to understand the underlying dependen- 
cies that define the shape of the IMF for galaxies and iden- 
tify which is the more fundamental driver of IMF variations. 
We have explored the variations with respect to SFRs, spe- 
cific SFRs and SFR surface density. These cases are de- 
scribed in !|5) !|6]2]and 

The best fit IMF slope (a) is determined for each SFR, 
specific SFR and SFR surface density sub-group of the 
three volume limited samples. This is essentially the IMF 
of the PEGASE model track closest to the region of highest 
data density. A library of model evolutionary tracks with 
input IMFs with high-mass slopes ranging from a = —1.45 
to —3.05 in increments of 0.1 is used to determine the best- 
fit IMF for each galaxy. The relationship between the IMF 
slope and mean property (i.e. SFR, specific SFR and SFR 
surface density) of objects in the respective sub-group is 



shown in Figure |13| The error in a is the standard devia- 
tion of the data in each respective sub-sample. The filled 
symbols in Fig ure|13| correspond to the bes t fit a slopes if 
Calzetti] (|2001[) /]Cardelli, Clayton fc Mathis| (|1989[) dust ob- 



scuration curves are used to correct the data, while the open 
symbols show the effect if the [Fischera fc Dopita| ( |2005[ ) 
dust obscuration curve is used instead. As shown in Fig- 
urepZ} dust obscuration corrections based on the IFischera] 
|fc Dopita| ( |2005[ ) curve seem to exaggerate the SFR-IMF 
trend. This can be seen explicitly in Figure 13. 

The trend to steeper IMF slopes with brighter Mr at 
a fixed SFR evident in Figure [T3{ ^a), is likely to be a conse- 
quence of higher luminosity systems having higher masses. 
For a fixed SFR these systems have smaller specific SFRs, 
and we have also demonstrated (Figures 1 1 1 1 and [Tsl^b) ) that 
galaxies with higher specific SFR have fiatter IMF slopes. 



Figure 13 b) highlights the result that, when accounting for 
the effect of galaxy mass, there is a surprisingly tight rela- 
tion between inferred IMF slope and specific SFR. It is pos- 
sible that there may also be metallicity effects contributing 
as well, although the population synthesis models suggest 
that these effects are likely to be small (Figure jsja). 

A least-squares fit to the SFR sub-groups of the high- 
est redshift volume limited sample gives 



Q ^ 0.361og(5Fi?) - 2.6. 



(13) 



Figure 13 c) shows the a versus SFR surface density 
(Esfr) relationship for the three volume limited samples. 
This relationship is as clean and, if anything, tighter than 
that seen with the specific SFR. This, again, is a conse- 
quence of accounting for galaxy size, this time through the 
proxy of surface area rather than mass, in refining the basic 
IMF-SFR relation. 

A least-squares fit to Esfa sub-samples of the highest 
volume limited sample gives 



Q^0.31og(EsFii)-1.7. 



(14) 



The models by Weidner fc Kroupa ( 2005 1 predict 



IGIMF slopes for galaxies assuming an underlying Salpeter 
IMF. The resultant IGIMF slope ranges for all galaxies 
from these models are always steeper than the slope of 
the underlying IMF, with low-mass galaxies {M^/Mq ^ 
10^) having a steeper and wider range of IGIMF slopes 
(—3.12 ^ ctiGiMF ^ —3.3) than high-mass galaxies 
(-3.07 > aiGiMF > -3.1 for M,/Mq « 10^°). The general 
trend predicted by the IGIMF models, steep IGIMF slopes 
for low-mass galaxies and vice versa, seems to be in agree- 
ment with the trend found here. However, the quantitative 
values for the slopes predicted by the IGIMF models for a 
Salpeter IMF are not, and in fact a much fiatter underlying 
IMF slope would be required to produce an IGIMF slope of 
around 2 (see also [Weidner, Kroupa fc Pflamm-Altenburg| 
|2010|. 

It is apparent that the trend to flatter IMF slopes is 
present as a function of each of SFR, specific SFR, and 
SFR surface density, although the differences between the 
volume-limited samples decrease most notably when SFR 
surface density is considered. It is reasonable, then, to infer 
that SFR surface density is most likely to be the underlying 
property on which IMF slope is primarily dependent, al- 
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log (SSFR) 



-2 -1.5 

log (T.SFRI 



Figure 13. The best fit IMF slope for each of the (a) SFR (Figure [ej, (b) specific SFR (Figure [TTjl and (c) SFR surface density 
(Figure T2]| su b-groups of the three volume limited samples used in this study. The filled symbols denote results when using the 
[Calzetti ( |2001[ l/ |CardeIli, Clayton &c Mathisj ( [l989| ) dust corrections and open symbols represent the [Fischera fc Dopitaj | |2005[ l dust 
corrections. The solid horizontal line indicates a Salpeter slope, the dot-dashed line indicates a Kroupa (2001) high-mass slope of 
a = —2.3 and the dashed line dictates the jBaldry fc Glazebrook] l |2003| IMF slope. 



though the trend with specific SFR is also very tightly con- 
strained between the different samples, and may be equally 
as significant. It is easier to imagine physical processes re- 
lated to SFR surface density (rather than to specific SFR) 
that could cause variation in the IMF. For this reason we 
propose here that SFR surface density (or more accurately 
the local space density of the SFR, quantified observation- 
ally as a surface density) is the underlying property gov- 
erning the solpe of the massive end of the IMF. 



[T995t also a gree with our conclusion. In addition, the mea- 
sured high SFRs (> lOOOM0yr~^) of submillimeter galajc- 
ies at 2 > 4 ([Michalowsk i, Watson fc Hjorth|20"l0 1 seem to 
require a much fiatter IMF slope of a « —1.5 to explain the 
data ( Baugh et al.|2005 l, which is again consistent with our 
conclusion. Even more intriguingly (although perhaps only 
coincidentally, in particular as SFR surface density is likely 
to be a more fundamental relation), the extrapolation of 
the linear relation between a and \og{SFR) measured here 
gives a « -1.52 for SFR = 1000. 



IMF OF THE MILKY WAY AND ITS 
NEIGHBOURS 



The current SFR of the Milky Way is 3 ± 1 Mgyr" ' ( Scalo 



1986| ). The Spitzer/IRAC GLIMPSE survey of the Galac- 
tic plane finds a total SFR of 0.68-1.45 Moyr"^ (iRobitaille 



fc Whitney||2010[ ) and using the data from WMAP and 
GLIMPSE surveys a Galactic SFR of 1.3 ± 0.2 Mgyr^^ 



is measured (Murray fc Rahman 20101. These SFRs are 



consistent with the Milky Way being placed close to the 
PEGASE model track with an input Salpeter IMF, lead- 
ing to the inference of a Salpeter-like IMF for the Milky 
Way. Although decades of observations of the stellar IMF 
within the Milky Way find slopes consistent with Salpeter, 
this is not inconsistent with our results. The fact that ex- 
ternal galaxies with SFRs similar to the Milky Way have 
similar IMF slopes is a valuable consistency check on our 
conclusions. 

The Milky Way does not have a constant star forma- 
tion history. The Milky Way SFR was higher in the past 
( Gilmore|[200l| ), in which case the early Milky Way would 
have had a flatter than Salpeter IMF (Figure 13 1. Studies 
of carbon-enhanced metal-poor stars report that the IMF 



of the early Milky Way was flatter than the present ( Tum- 



linson 2007 Lucatello et al. 2005), again consistent with 



our results. 

The measured low SFRs of the Magellanic Clouds 
(O.14M0yr~^ for the Large Magellanic Cloud and 
0.015 M0yr~'^ for the Small Magellanic Cloud as given by 
Murray fc Rahman|2010 \ and the measured IMF slopes for 
their stellar clusters being Salpeter or steeper ( Massey et al. 



9 SALPETER IMF VS OTHER WIDELY 
ADOPTED IMFS 

The analyses presented in this study use an IMF with a 
Salpeter high-mass slope. However, there are other common 
forms of IMF in the literature. Here we explore how the 
model evolutionary paths vary if an IMF with a different 
analytical form or slope is assumed. As Figure[T4] shows, 
the commonly used IMFs in the literature produce model 
evolutionary paths that lie below the Salpeter IMF model 
track used in this study. Hence, unless the high-mass slope 
is adjusted appropriately, none of the common IMFs can 
produce an evolutionary track that describes the high SFR 
galaxies in the GAMA sample. 



10 THE DEGENERACY WITH RESPECT TO 
TURNOVER MASS 

Here we explore the degeneracy with respect to turnover 
mass. The turnover mass represents the mass at which the 
two-part power-law IMF turns over. 



1 |Kroupa| |200l] |A IMF: a = -1.3 for 0.1 < M/Mq < 0.5 and 
a = -2.3 for 0.5 < M/Mq < 120. 

|Kroupa| l |200ll )B IMF: a = -1.8 for 0.1 < M/Mq < 0.5 and 
a = -2.7 for 0.5 < M/M@ < 1 and o = -2.3 for 1 < M/Mq < 
120. 
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Figure 15. (a) The effect of varying the turnover mass of the IMF in the case of two-part power law IMF and (b) three-part power 
law IMF. See the text for further details. 




-2.35, 1 < M/Mq < 120. 



(19) 



Figure 14. Solid lines from top-to-bottom: [Kroupa] | |200l| A^ 
and |Kroupa| |2C)01| B^ IMFs produce model evolutionary tracks 
that overlap with one another. The rest of the evolutionary tracks 
are generated by using Scalo| (|1998^, |Kennicutt| (|1983^ , [Kroupa,] 
[Tout fc Gilmore| | [T993 I and [Miller fc Scalo| | |1979| l IMFs. The 
dashed lines represent the three evolutionary model tracks from 
previous figures. 



The solid lines (from top-to-bottom track) in Fig- 



ure 15 a) are: 
a= -1.3, 0.1 



a = 



a = 



-1.3, 
-2.35 
a= -1.3 



a = 



< M/Mq < 10 
2.35, 10 < M/Mq < 120. 
-1.3, 0.1 < M/Mq < 5 
-2.35, 5 < M/Mq < 120. 
0.1 < M/Mq < 3 

3 < M/Mq < 120. 
0.1 < M/Mq < 2 
-2.35, 2 < M/Mq < 120. 
-1.3, 0.1 < M/Mq < 1 



(15) 
(16) 
(17) 
(18) 



The degeneracy with respect to turnover mass in the 
IMF is such that our results could be explained by invoking 
not a change in the high-mass slope of the IMF, but a 
progressive increase in the stellar mass at which the IMF 
turns over at the low-mass end. As the mass of the IMF 
turnover increases the model tracks shift upwards in the 
same manner as seen when making the high-mass slope 
flatter. This is illustrated in Figure 15 a), which shows that 



extending the low mass range up to 10 Mq produces a 
model track that is similar to the top dashed track. 



Figure 15 b) presents the variations of the model tracks 
for different combinations of slopes and mass ranges of 3- 
part power law IMFs. 

The solid lines (from top-to-bottom track) in Fig- 
ure 15 b) are: 



= -1.3, 



0.1 < M/Mq < 0.5, 



-1.6, 0.5 < M/Mq < 10 and 



10 < M/Mq < 120. 
-1.3, 0.1 < M/Mq < 0.5, 
-2.0, 0.5 < M/Mq < 5 and 
5 < M/Mq < 120. 
, 0.1 < M/Mq < 0.5, 
-2.0, 0.5 < M/Mq < 10 and 
2.35, 10 < M/Mq < 120. 



-2.35, 



-2.35 
-1.3 



(20) 



(21) 



(22) 



While it is possible that our results could be explained 
by a modification of the IMF such that the turnover mass 
increases with SFR, rather than our claimed flattening of 
the high-mass slope, this seems unlikely. To reproduce the 
high-SFR GAMA systems would require a turnover mass of 
~ 10 M0 , which seems surprisingly high. This could still be 
a possibility, however, and is included here as an alternative 
explanation for completeness. 
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11 SUMMARY 

We have used ~ 33 000 galaxies from the GAMA survey 
to confirm that the IMF does not appear to be universal. 
We have shown that the stellar IMF within galaxies has a 
strong variation with galaxy SFR. 

This result is consistent with many recent studies that 
suggest an evolving or varying IMF as a solution to the 



observed discrepancies. Many authors ( Hopkins & Beacom 



20061 


Wilkins, Trentham & Hopkins 2008b| Wilkins et al. 


2008a 


Fardal et al 


2007 Perez-Gonzalez et al. 2008 van 


Dokkum|2008 


Dave 


20081) have suggested an evolving IMF 



in order to reduce the discrepancy between the observed 
stellar mass density of the Universe and that implied by the 
cosmic star formation history. Accordin g to [Wilkins, Tren-| 



tham & Hopkins|(|2008b'); 'Wilkins et al.|(|2008al) a "cosmic" 



IMF with a = -2.15 ( ;Baldry fc Glazebrook||2003[ ) would 



solve the discrepancy between the two quantities at low red- 
shift {z < 0.7) but an IMF that is still flatter is required 
for z > 0.7. Given that high redshift sources also tend to 
be the high SFR objects, our results predict the IMF of 
high-redshift galaxies to be very fiat. Additionally, studies 
looking at the evolution of mass-to-light ratios and colours 
of galaxies (e.g. [van Dokkum|[2008| show that the models 
that best fit the observations are those that assume a flatter 
IMF at high redshift. Galactic chemical models of |Calura fc| 
|Menci| p009) require a flatter IMF for massive galaxies in 
order to correctly predict their observed metallicities. We 
can now provide an explanation for these suggestions, with 
the finding that the IMF has a strong underlying depen- 
dence on the host galaxy SFR. 
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